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Abstract 

The  practical  application  of  the  quantum  1/f  effect  to  quartz  resonators  and  to  infrared  detectors 
present  in  this  report  allows  us  for  the  first  time  to  understand  and  to  extend  the  stability  limits  of 
quartz  resonators.  It  also  explains  1/f  noise  in  most  semiconductor  devices  and  in  infrared  detectors 
considered  in  this  report  in  the  presence  of  radiation,  aithough  1/f  noise  has  been  already  successfully 
pushed  below  background  in  PtSi  on  p-type  Si  Schottky  diodes  at  RADC/Hanscomb.  A  fundamental 
breakthrough  was  performed  through  the  first  direct  derivation  of  the  coherent  quantum  1/f  effect 
from  a  special  quantum-eiectrodynamic  propagator,  and  from  the  author's  general  sufficient  1/f  chaos 
criterion  presented  in  the  previous  yearly  report.  Finally,  the  quantum  1/f  cross-correlations 
derived  by  the  author  have  been  used  to  recalculate  and  to  graph  the  quantum  1/f  mobility  fluctuations 
in  Si  and  GaAs  samples  as  a  function  of  temperature  and  doping,  in  good  agreement  with  the 
measurements  of  Tacano  in  Japan  and  Hooge  in  the  Netherlands. 
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1.  INTRODUCTION 


The  progress  achieved  this  year  includes  the  derivation  of  the  Coherent  Quantum  1/f  Effect  from 
a  special  cpiantum-electrodynamicai  propagator  known  since  1975,  and  from  the  sufficient  1/f 
Oriterion  bitroduced  in  the  previous  Annual  Report  (Sec.  11  and  III.  respectively).  It  also  includes  the 
successful  practical  application  of  the  quantum  1/f  theory  to  quartz  resonators  and  the  verification  of 
the  results  at  NiST-Boukter  with  the  cooperation  of  Fred  Walls  (Sec.  IV).  it  further  includes  in  Sec.  V 
the  practical  s^ication  of  the  quantum  1/f  theory  to  infrared  detectors  in  the  presence  of  radiation,  in 
particuiar.  a  proof  of  the  absence  of  quantum  1/f  noise  in  the  process  of  carrier  photogeneration  in 
semiconductors,  it  finally  includes  in  Sec.  Vi  the  recalculation  and  first  time  graphical  representation 
of  cpiantum  1/f  mobility  fluctuations  in  Si  and  GaAs  samples,  based  on  the  author's  new  cross- 
correlation  formula,  as  a  function  of  temperature  and  doping  level. 

We  also  mention  the  NSF-sponsored  organization  by  the  author  of  the  "V.  international  van  der 
Ziei  Symposium  on  Quantum  1/f  Noise  and  Other  Low-Frequency  Fluctuations”  at  the  University  of 
Missouri-St.  Louis  on  May  22-23,  1992  with  a  100%  larger  participation  than  ever  before,  and  the 
ongoing  organization  by  the  author  of  the  *XII  Int  Conf.  on  Noise  in  Physical  Systems  and  1/f  Noise”  in 
St.  Lois,  for  Aug.  16-20,  1993.  Finally,  we  mention  the  aeation  of  an  institute  of  Molecular 
Electronics  at  the  University  of  Missouri  St.  Louis  with  the  author's  participation  in  September  1991. 
My  ooiiaboralors  have  been  F.  Wails,  E.  Bernard!,  T.  Chung,  A.  FOrst-Chung,  Y.  Zhang,  Xuewei  Hu,  Jian 
Xu,  and  I.  ProleHto.  These  results  will  be  briefly  presented  below. 

II.  COHERENT  QUANTUM  1/f  CHAOS 

Conventional  quantum  1/f  fluctuations  of  physical  cross  sections  and  process  rates  have 
been  introrkjced  by  us  as  a  fundamental  infrared  divergence  phenomenon  in  1975  [1].  Some  of 
the  subsequent  publications  I2]-[141  have  shown  this  new  effect  to  be  unaffected  by  the 
presence  of  the  thermal  rarfiation  background  [4],  [5],  some  have  derived  it  with  wave  packets 
(11],  including  a  finite  mean  free  path  [9],  in  second  quantization  [13],  with  the  Keidysh- 
Schwkiger  method  [12],  or  in  the  Van  Hove  weak  interaction  limit  [14],  some  derived  its 
chan^teristic  functional  [7],  or  applied  it  to  the  calculation  [10]  of  mobility  and 
recombination  speed  fluctuations  In  semiconductors  and  semiconductor  devices.  Others  [15]  - 
[27]  verified  the  new  effect  experimentally  and  successfully  applied  it  to  electronic  devices. 

The  present  paper  derives  a  related  fundamental  effect  which  we  call  the  coherent  1/f 
effect,  with  eMmentary  methods  of  quantum  electrodynamics  and  non-relativistic  many-body 
ttieory.  Our  derivation  uses  the  new  picture  introchiced  by  Doilard  [28]  and  generalized  by 
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Kuiish  and  Fadde’ev  [29]  and  later  by  Zwanziger  [30],  in  agreement  with  earlier  work  by 
Chung  [31]  and  by  KIsbie  [32].  in  this  new  picture,  the  asymptotic  Coutomb  inter£K:tion  is 
included  in  the  unperturbed  Hamiltonian  rather  than  in  the  perturbation  part.  This  ieads  to  a 
more  ocmptBx  physical  free  particle  notion  which  includes  a  coherent  photon  doud,  and  replaces 
the  pole  in  the  propagator  with  a  branch  point,  it  aiso  leads  to  a  smeared-out  mass  shell.  Using 
this  picture,  we  can  neglect  the  remaining  part  of  the  interaction  if  we  limit  ourselves  to  the 
asymptotic  region  of  large  distances  and  times,  which  are  important  in  the  case  of  1/f  noise. 

For  N  electrons  in  a  Fermi  sphere  shifted  in  momentum  space  by  a  vector  po  and 
occupying  N/2  orbitals  the  propagator  derived  by  these  authors  [30]  can  be  reduced  for 
large  time  components  of  x'-x  to  the  non-relativistic  form 

-i«&o|Tv8*(x')¥l‘(x)|«»o>  ■  8ss‘  Gs(x*-x) 

-  (i/V)I{expi[p(r-r*)-p2(t.f)/2ml/hr}np.s 

P 

x{-lp(r-r')/hr+i(m2c2+p2)i/2(t.f)(c/K)}“/*.  ( 1  ) 

Here  a«e2/llc»1/137  is  Sommerfeld's  fine  structure  constant,  np,s  the  number  of  electrons  in 
the  state  of  momentum  p  and  spin  s,  m  the  rest  mass  of  the  fermions,  Sss*  the  Kronecker 
symbol,  c  the  speed  of  light,  x«(r,t)  any  space-time  point  and  V  the  volume  of  a  normalization 
box.  T  is  the  timeK>rdering  operator  which  orders  the  operators  in  the  order  of  decreasing 
times  from  left  to  right  and  multiplies  the  result  by  (-1)^,  where  P  is  the  parity  of  the 
permutation  required  to  achieve  this  order.  For  equal  times,  T  normal-orders  the  operators, 
i.e.,  for  t-f  the  left-hand  side  of  Eq.  (1)  is  i<Oo|y|(x)vs'(x')|Oo>'  The  state  Oq  of  the  N 
electrons  is  described  by  a  Slater  determinant  of  single-particie  orbitals. 

Consider  first  the  case  M'  for  simplicity,  although  only  the  case  of  large  t-f  can  be 
expected  to  be  experimentaPy  applicable.  The  pair  correlation  function  can  then  be  decomposed 
as  follows 

<^ol¥t(x)vl*(x’)V8*(x’)y»(x)|<Po>  -  <»PoIvI(x)v*(x)|«Po><»&oIvI*(x')Y8’(x  )|«o> 

*  <‘frolYt(x)V8'(x‘)lOo>«PolVl'(x')v8(x)|<Po>  •  (  2 ) 

The  first  term  can  be  expressed  in  terms  of  the  partide  density  of  spin  s,  n/2  -  N/2V  - 
<^oiYt(x)Y»(x)l^o>>  whHe  the  second  term  can  be  expressed  in  terms  of  the  Green  function 
(1)  in  the  form 

Am*(x-x*)  ■<®olYl(x)Yt*(x')Y8’(x’)Y8(x)|«o>  -(n/2)2  Q8(x’-x)G8(x-x’).  (3) 
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The  'relative*  autocorrelatton  hjnction  A(x-x')  describing  the  normalized  pair  oorreiation 
independent  of  qpin  te  obtsrined  by  cfivkftig  by  and  summing  over  s  and  s' 

A(x-x*)  -  1  +  (1/n2)5:Gs(x-x')G8(x'-x) 

8 

-  1  -  (1/N2)Z  I{expli(p-p')(r-r')/hri}np.8np-.s 

s  pp' 

x|p(r-r')/hr|«'*|p'(r-r')/hr|“^*  (4) 

Here  we  have  used  Eq.  (1).  The  iow-wavenumber  part  A|  of  this  relative  density 
autocorrelation  function  is  given  by  the  terms  with  p-p’. 

Ai(x-x')  -  1  -  (1/n2)Z  Inp.8lP(r-r')/h'|2“/*  (5) 

8  P 

-  1  -  (2/N2)[V/(2kK)3]  Jd3p|(p  +  po)(r-r’)/K’|2«/« 

P<Pf 

a  a  a 

-  1  -  (2/N2){V/(2xK)3I  /(*>,  !dp2  idp3l(P3+Po3)(r-r')/K|2“'« 

-a  -a  -a 

-  1  -  (2/N2)[V/(2jiK)3]8a3la(r-r')/h'|2o/« 


1  -  (1/N)[(ii/6)i/3pp|r.r'|/lf]2o/K. 

for  Pf»Po3! 

(6) 

1  -  (1/N)|po(r-r’)/h|2«/" 

for  Pf«Po3! 

(6*) 

In  view  of  the  smallness  of  ZoJk,  for  Pf»Po3  we  have  integrated  in  cartesian  coordinates, 
approxbnating  the  Fermi  sphere  by  a  cube  of  side  2a  with  a-  (ic/6)^^Pf-  The  result  is 
fxracticaliy  independent  of  po  and  of  Po3-Po('’*i’')/|>^-f'l-  The  factor  (jc/6)2o/3k  can  be 
negfocted.  For  Pf«Po3  we  used  the  mean  value  theorem  for  estimating  the  integral  over  d3p  in 
spherical  coordinates,  in  ail  cases  the  autocorrelation  decreases  very  slowly  from  1  when  |r- 
r*!  is  focreased  to  very  large  vsdues.  Writing  the  rectangular  bracket  in  Eq.  (6)  as  an 
exponentiai  funcdon  of  its  logarithm,  expandina  the  resulting  exponential,  and  keeping  only  the 
^  term,  we  obtafo  with  pp/Kakp 

Ai(x-x')  -  {1-2/N)  ♦  (6/*)2«/3i»{ii/kF|r-r'|l2»/*}/N 
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-  (1-2/N)  +  (6/ii)2«/3"[1/kFl2«'*(2o/Nw)Jcos[k|r-r'|ldk/ki-2«^*. 

oo 

-{N-2  +(2a/x)J[k/kF]2a^*cos[k|r-r*|ldk/k}/N  for  Pf»Po3:  (7) 

OO 

-{N-2  +(2a/ic)J[k/kol2«'*cos[kpo(r-r*)/Pol  clk/k}/N  for  pF«po3;  (7*) 

Here  we  have  used  a  well-known  Fourier  integral  [35]  and  we  have  introduced  ko9po3/li. 
According  to  the  Wiener-Khintchine  theorem,  the  coefficient  of  the  cos  gives  the  spectral 
density.  To  get  it  for  the  fractional  fluctuations  Sn/n,  we  divide  by  the  constant  term  N-2 

S8n/n(k)  -  [2a/iik(N-2)][k/Kl2«/*.  (8) 

where  K-kf  for  Pf»Po3  as  in  the  case  of  metals  with  spherical  wave  symmetry,  and  K-ko  for 
Pf«Po3<  Although  inapplicable,  this  pure  1/k  spectrum  is  the  wave-number  equivalent  of  the 
coherent  quantum  1/f  noise  derived  earlier  [33],  [34]  in  excellent  agreement  with  the 
experiments  on  large  electronic  devices  [15],  [25].  Due  to  2aA«1  the  second  factor  is 
practically  unity  and  of  no  importance,  except  for  eliminating  the  logarithmic  divergence  from 
the  spectral  integral.  This  wave  number  spectrum  also  entails  a  1/f  frequency  spectrum 
obtained  by  writing  (fi</k-df/f,  as  was  shown  in  detail  in  a  previous  paper  [13].  For  equal 
times  our  result  can  not  be  expected  to  be  valid,  due  to  the  asymptotic  character  of  Eq.  (1).  We 
shall  now  derive  the  1/f  spectrum  directly  below, 
if  Eq.  (2)  is  replaced  by 

<‘»o|Tirt(x)v8(x)vJ‘(x’)Y8’(x’)|<6o>  -  «»olvl(x)y8(x)!<&o>«^olvt’(x’)v8’(x*)|Oo> 

-  «»o|Tv8’(x')y|(x)lOo><^o|Ty8(x)yi'(x’)l^o>.  (9) 

Eq.  (3)  remains  the  same,  except  for  the  middle  part  which  is  replaced  by  the  left  hand  side  of 
Eq.  (9).  Eq.  (4)  becomes  now 

A(x-x*)  -  1  -  (1/n2)2:Q8(x-x*)G8(x’-x) 


1  *  (1/n2]I  I{expl[(p-p-)(r-r*)-(p2.p‘2)(t.f)/2m]/h'}np.8np*.8 

•  PP' 
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x{p(r-r’)/h'-(m2c2+p2)i/2(t.f)(c/tf)}«/* 

X  {p  •(r-r’KhT-  ( m2c2+  p  -2)  i  /2(t.f)(c/tf)}“^*. 


(10) 


H«r«  we  have  used  agiAi  B:|.  (1).  The  low-frequency  and  low-wavenumber  part  A|  of  this 
relative  density  autocorrelation  function  is  aiso  given  by  the  terms  with  p-p*. 


Ai(x-x')  -  1  -  (1/n2)5:  Snp.s 

s  p 

x{p(r-r’)/h'-(m2c2+p2)i/2(t.f)(c/K)}2«/"  (11) 

-1  -(2/N2)lV/(2xK)«  Jd»p|(p  +  po)(r-f)/tf-[m2c2+(p  +  po)2]’'2(t-t’)(c/R)|2a/« 

fxpF 

-  1  -  (1/N)|po(r-r‘)/h'’-mc2x/K|2a/*  for  pF«|po3-mc2T/z|.  (12) 


Here  we  have  used  the  mean  value  theorem,  considering  the  lain  power  as  a  siowiy  varying 
function  of  p  and  neglecting  Po  in  the  coeffident  of  t  ■t-t',  with  ZB|r-r'|.  Writing  the  power 
again  as  an  exponential  function  of  its  iogarithm.  expanding  the  resuiting  exponentiai,  and 
keeping  only  the  first  term,  we  obtain  with  (tf/mc2)|po(r-r’)/h’-mc2t/K|Me«|T-po(r* 
r')/mc2j 


Ai(x-x')  -  1-[(mc2/lf)0]2o/K/N 


(1-2/N)  +  [lf/mc2]2o/x(2o/Nx)Jcosla)e]da)/<Di 


{N-2  +  (2o/K)J^[K<o/mc2]2«^*COslwej  d«()/a)}/N. 


-20/ ji 


(13) 


According  to  ttte  Wiener-Khintchine  theorem,  the  coefficient  of  the  cos  gives  the  spectrai 
density.  To  get  it  for  the  fractional  fluctuations  Sn/n,  we  divide  by  the  constant  term  N-2 

Sin/nik)  -  pa/ic®(N.2)](K©/mc2l2o/K.  (14) 


This  rMuit  is  also  applicable  for  the  particle  current  fluctuation  spectrum,  indeed,  for  current 
dwishy  fhictuadons  fij  we  include  a  (K/mi)V  in  front  of  each  of  the  two  y  operators  in  Eq.  (9), 
a  factor  ppVpo^  in  Eq.  (10)  after  the  summation  signs,  a  factor  (p/po)^  in  the  first  form  of  Eq. 
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(t1),  a  fector  (p-t'Po)^/Po^  in  tha  second  form,  and  no  changes  in  Eqs.  (12)-(13).  Eq.  (14) 
becomes 

Sgjflfkl  -  f2a/rcmfN.21UKQ)/mc2l2ayic.  (15) 

This  result  coincides  with  our  earlier  theoretical  result  for  coherent  quantum  1/f  noise  if  we 
replace  N  with  N-2.  The  validity  of  this  equation  is  restricted  to  low  frequendes  and  wave- 
numbers.  This  equation  is  in  excellent  agreement  with  mobility  and  diffusion  1/f  noise  in  large 
devices. 

Rnally,  we  consider  the  errors  caused  by  ttie  neglect  of  higher  order  terms  in  the 
expansion  of  the  exponential  functions  resulting  from  Eqs.  (6)-(6')  and  (12).  For  a  thermal 
electron  and  r«1cm  in  Eq.  (6')  we  get  kr-IO^^  and  (2a/ic)ln(kr)-0.12,  yielding  an  error  of 
12%.  For  t-IO^^s  in  Eq.  (12),  which  is  the  age  of  the  universe,  we  get  an  error  of  40%  as  the 
upper  limit.  Therefore  a  more  exact  treatment  is  of  some  interest.  Using  the  identity  [35] 

oo 

^lai*  *  ['(2a/ic)  /co*2®^*COS(0(D)d(io/<fll 
«»0 

x{cosa  +  (2o/x)]£(0(a«)2n-2a/«[(2n)!(2n.2a/J^)]•■'}•^  (16) 

n>0 


with  arbitrarily  small  cutoff  <i»o>  we  obtain  from  Eq.  (12)  the  exact  form 

A|(x-x')  «  1  +  [(2o/xN)  J(mc2/Ka))2a/*cos(0a))d<ii)/a>l 

•• 

x{coso  +  (2a/jt)5^(0<Bo)2n-2o/K[(2n)l(2n-2a/it)]*1}*1.  (17) 

n-0 

This  would  indcate  a  <0*i'2a/K  spectrum  and  a  1/N  dependence  of  the  spectrum  of  fractional  n 
and  j  fluctuafions,  if  we  neglect  the  curly  bracket  in  the  denominator  which  is  close  to  unity  for 
very  smaH  <oo-  We  thus  realize  that  the  unusuai  N-2  dependence  in  Eqs.  (14)-(15)  is  caused 
by  the  forced  fotroduction  of  the  integrate  spectmm  in  place  of  the  spectrum. 

Due  to  the  smdiness  of  a  both  forms  coincide  in  practical  applications.  Eq.  (15)  for  the 
coherent  QED  chaos  process  in  electric  currents  can  thus  be  written  also  in  the  form 

Sk|/i(k)  -  f2a/tccoNlfmc2/K<fli2o^-2a/KOON  -  0.0046S/CQN.  (18) 
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This  result  derived  directly  earlier  [33],  [34],  is  in  excellent  agreement  with  the 
measurements  [15],  [25],  in  large  {see  [34]  for  a  definition  of  "large”  or  "extended”,  and  for 
an  interpolifiion  with  the  conventional  quantum  1/f  chaos  effect}  devices  such  as  large  n+p  Hgi . 
xCdxTe  infrared  detector  diodes.  It  is  also  dose  to  the  empirical  value  of  0.002/(dN  observed 
earlier  by  Hooge  [36]  in  semiconductors  and  metals.  Being  observed  in  the  presence  of  a 
constant  applied  field,  these  fundamental  quantum  current  fluctuations  are  usually  interpreted 
as  mobility  fluctuations. 

ill.  APPUCATION  OF  THE  SUFFICIENT  1/f  CRITERION  TO  QUANTUM  1/f  CHAOS 


The  nonlinearity  causing  the  1/f  spectrum  of  turbulence  in  both  semiconductors  and 
metals  is  caused  by  the  reaction  of  the  field  generated  by  charged  particles  and  their  currents 
back  on  themselves.  The  same  nonlinearity  is  present  in  quantum  electrodynamics  (QED), 
where  it  causes  the  infrared  divergence,  the  infrared  radiative  corrections  for  cross  sections 
and  process  rates,  and  the  quantum  1/f  effect.  We  shall  prove  this  on  the  basis  of  the  suffident 
criterion  for  1/f  spectral  density  in  chaotic  systems,  derived  in  the  previous  annual  report. 

Consider  a  beam  of  charged  particles  propagating  in  a  well-defined  direction  which  we 
shall  call  the  x  direction,  so  that  the  one-dimensionai  SchrOdinger  equation  describes  the 
longitudinal  fluduations  In  the  concentration  of  particles.  Considering  the  non-relativistic  case 
which  is  encountered  in  most  quantum  1/f  noise  applications,  we  write  in  second  quantization 
the  equation  of  motion  for  the  Heisenberg  field  operators  v  of  the  in  the  form 

Ih8v/3t  -  (1/2m)[-ihV  -(e/c)  A]2y,  (19) 


With  the  non-relativistic  form  J  «  •ihv*Vy/m  +  hermitic  conjugate,  and  with 


A(x,y,2,t) 


(20) 


we  obtain 

ihdv/dt  -  (1/2m)[-inv  -(eri/2c2mi) 


(21) 


At  very  low  freqi^ndes  or  wavenumbers  the  last  term  in  rectangular  brackets  is  dominant  on 


the  r.h.8.,  leacfing  to 
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Ihdv/dt  -  (-1/2m)[(eri/2c2m)  ^dxj%.  (22) 

For  X  replacod  Xx,  and  x*  replaced  by  Xx',  we  obtain 

ihdy/dt  «  (-1/2m)[(eri/2c2m)J^'’^  ^X2dx']%  »  X^Hy  *  X“PH\|r.  (23) 

This  satisfies  our  homogeneity  criterion  with  p— 2.  Our  sufficient  criterion  only  requires 
homogeneity,  with  any  value  of  the  weight  p,  for  the  existence  of  a  1/f  spectrum  in  chaos. 
Therefore,  we  expect  a  1/f  spectrum  of  quantum  current-fluctuations,  i.e.,  of  cross  sections 
and  process  rates  in  physics,  as  derived  in  detail  in  Sec.  11  above.  This  is  in  agreement  with  the 
well-known,  and  experimentally  verified,  results  of  the  Quantum  1/f  Theory. 

in  conclusion,  we  realize  that,  both  in  classical  and  quantum  mechanical  nonlinear 
systems,  the  limiting  behavior  at  low  wave  numbers  is  usually  expressed  by  homogeneous 
functional  dependences,  leading  to  fundamental  1/f  spectra  on  the  basis  of  our  criterion. 


IV.  QUANTUM  1/f  FLXJKER  OF  FREQUENCY  IN  QUARTZ  RESQNATQRS:THEQRY  AND  EXPERIMENT 

IV.  1.  Introduction 

Flicker  of  frequency  noise  is  an  important  characteristic  of  quartz  resonators  which  limits 
their  stability  and  determines  their  utility  for  most  applications.  The  1/f  contribution  to  frequency 
stability  is  best  obtained  by  observing  the  stability  over  a  range  of  measurement  times  in  the  time 
domain,  or  over  a  range  of  frequencies  in  the  frequency  domain.  From  the  extended  data  one  can  fit  a 
flicker  of  frequency  model  to  the  data  that  excludes  the  contributions  from  random  walk  frequency 
modulation  and  from  the  drift  present  in  many  resonators  and  oscillators  (Fig.  1). 

One  of  the  first  who  systematically  studied  the  1/f  noise  as  a  function  of  geometry,  temperature 
and  Q-factors  was  Gagnepain  [37].  He  noticed  empirically  that  the  1/f  part  of  the  spectral  density  of 
fractionai  frequency  fluctuations,  Sy(f).  varied  as  Q'^  for  resonators  between  1  and  25  MHz.  As  the 
temperature  changes,  the  Q-^ctor  of  a  resonator  changes.  This  allows  us  to  exclude  the  effect  of  many 
other  tabors.  Additional  work  by  Parker,  however,  showed  that  the  data  from  both  bulk  acoustic  wave 
(BAW)  and  surface  acoustic  wave  (SAW)  devices  could  be  roughly  fit  to  the  the  same  model  [38],  if  one 
assunras  a  dependence  for  the  phase  noise  kisted  of  Sy.  The  fit  is  shown  on  Fig.  2. 

From  a  theoretical  point  of  view,  fundamental  work  by  this  author  [39],  inspired  by  the 
quanttjm  1/f  theory  [40],  has  derived  the  Q‘^  depenence  of  Sy  for  all  resonators  from  first  principles 


afre«ly  1979.  Work  on  many  systems  other  than  quartz  has  yielded  very  good  quantitative  agreement 
between  theory  and  experimental  data  for  quantum  1/f  noise  [41].  This  paper  refines  the  previous 
theoretical  work  on  1/f  frequency  noise  in  quartz  to  suggest  a  better  framework  for  predicting  the 
level  of  1/f  frequency  noise  in  quartz  resonators  over  a  wide  range  of  frequencies  and  Q-factors. 

IV.  2.  Theory  of  1/f  Frequency  Noise  in  Quartz  Resonators 

According  to  the  general  quantum  1/f  formula  [40],  r-2Sr(f)*2aA/f  with  a-e^/fic  -1/137 

and  A-2(AJ/ec)^/3n  is  the  quantum  1/f  effect  in  any  physical  process  rate  r.  Setting  J-dP/dt-P 
where  P  is  the  vector  of  the  dipole  moment  of  the  quartz  crystal,  we  obtain  for  the  fluctuations  in  the 
rate  r  of  phonon  removal  from  the  main  resonator  oscillation  mode  (by  scattering  n  a  phonon  from  any 
other  mode  of  average  frequency  <(■)>)  of  the  crystal,  (or  via  a  two-phonon-process  at  a  crystail  defect 
or  impurity,  involving  a  phonon  of  average  frequency  <(d‘>)  the  spectral  density 

Sr(f)  -  r24o(AP)2/3jte2c2.  (24) 

where  (AP)2  is  the  square  of  the  dipole  moment  rate  change  associated  with  the  process  causing  the 
removal  of  a  phonon  from  the  main  oscillator  mode.  To  calculate  it,  we  write  the  energy  W  of  the 
interacting  resonator  mode  <(a>  in  the  form 

W-nri<a)>  -  2(Nm/2)(dx/dt)2.(Nm/e2)(e  dx/dt)2-(m/Ne2)e2(P)2;  (25) 

The  factor  two  includes  the  potential  energy  contribution.  Here  m  is  the  reduced  mass  of  the 
elementary  oscillating  dipoles,  e  their  charge,  e  a  polarization  costant,  and  N  their  number  in  the 
quartz  crystal,  /^plying  a  variation  An-1  we  get 

An/n  -  2|AP|/1P|,  or  AP-P/2n.  (26) 

Solving  Eq.  (25)  for  P  arxl  substituting,  we  obtain 

|AP|  -  (Nh<a»/n)1/2(e/2£)  (27) 

Substituting  aP  into  Eq.  (24),  we  get 

r-2Sr(f)  -  Nan<a»/3nKmc2fe2  -  A/f .  (28) 
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This  result  is  applicable  to  the  fluctuations  in  the  toss  rate  r  of  the  quartz. 

The  corresponding  resonance  frequency  fluctuations  of  the  quartz  resonator  are  given  by^ 

w"^Sa)(f)  -  (1/4Q)(A/f)  -  Nah<(o>/12n«mc2fe2Q4,  (29) 

where  Q  is  the  quality  factor  of  the  single-mode  quartz  resonator  considered,  and  <a)>  is  not  the 
circular  frequency  of  the  main  resonator  mode.  (Oq.  but  rather  the  practically  constant  frequency  of  the 
average  interacting  phonon,  considering  both  three-phonon  and  two-phonon  processes.  The 
corresponding  AP  in  the  main  resonai^^r  mode  has  to  be  also  included  in  principle,  but  is  negligible 
because  of  the  very  large  number  of  phonons  present  in  the  main  resonator  mode. 

Eq.  (29)  can  be  written  in  the  form 

S(f)  -  pV/fQ4.  (30) 

where,  with  an  intermediary  value  <a)>3:  10®/s,  with  n=kT/h<a)>,  T=300K  and  kT=4  10^^, 

p  »  (N/V)ah«o>/12n7ce2mc2  =  1022(1/137)(10-271  o8)2/i2kTji1  0*27  9  1020  =  1.  (31) 

The  form  of  Eq.  (30)  shows  that  the  level  of  1/f  frequency  noise  depends  not  only  as  Q*^  as 
previously  proposed,  but  also  on  the  oscillation  frequency  or  the  volume  of  the  active  region.  This 
model  qualitatively  fts  the  data  of  Gagnepain  [37,48  ]  where  he  varied  the  Q-factor  with  temperature 
in  the  same  resonator  (but  not  frequency  or  volume). 

The  model  also  provides  the  basis  for  predicting  how  to  improve  the  1/f  level  of  resonators, 
beyond  just  improving  the  Q-factor,  which  has  been  known  for  many  years.  Since  the  level  depends  on 
active  volume,  one  should  use  the  lowest  overtone  and  smallest  diameter  consistent  with  other  circuit 
parameters. 

IV.  3.  Experimental  Measurements  and  Analysis  of  1/f  Noise  in  Quartz  Resonators 

The  level  of  1/f  frequency  noise  in  quartz  resonators  has  been  measured  using  phase  bridges  and 
complete  oscillators  [37,38,42-47].  Unfortunately  much  of  the  data  in  the  literature  is  unusable  for 
mocteling  because  the  unloaded  Q-factor  is  unknown.  (Qur  case  is  even  more  restrictive  because  we  also 
need  to  know  the  electrode  size).  The  phase  bridge  approach  has  the  advantage  that  the  unloaded  Q- 
factor  can  be  easily  measured  and  the  noise  in  the  measurement  electronics  can  be  evaluated 
independent  of  the  resonator.  If  resonator  pairs  are  used,  the  noise  of  the  driving  source  can  generally 
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te  no(^o^ed  and  the  pair  can  operate  at  virtuaity  any  frequency  [45].  The  oscillator  approach  makes  it 
possible  to  compare  many  different  resonators  one  at  a  time.  The  noise  of  imfividual  oscillators  can  be 
derived  by  measuring  the  phase  noise  between  3  oscillators  [48]. 

Figure  3  taken  from  [38]  is  one  of  several  studies  showing  that  the  1/f  level  is  virtually 
independent  of  the  loaded  Q-factor.  This  is  in  complete  agreement  with  the  theoretical  model.  In 
practicai  oscillators  there  is  a  dependence  on  loaded  Q-factor  only  when  the  phase  notee  of  the 
sustainkig  electronics  contributes  to  the  overall  noise  level. 

We  have  analyzed  l/f  frequency  noise  as  a  function  of  unloaded  Q,  volume  under  the  electrodes, 
and  frequency.  For  a  given  resonator  geometry  and  manufacturer,  we  have  taken  the  best  values  of 
Sy(Q  reported  in  an  attempt  to  remove  the  effects  of  poor  crystals  or  electronics.  In  Fig.  4  we  have 
taken  all  of  the  precision  data  available  with  unloaded  Q-factor,  electrode  volume,  and  frequency 
stability  and  plotted  them  according  to  the  three  models.  Except  for  the  2.5  MHz  resonator  where  Qvq  > 
0.95x10^3,  the  Qvo  product  for  all  resonators  is.  near  1.2x1013  (this  is  close  to  the  material  limit  for 
AT  and  SC  cut  resonators).  The  curve  labeled  Ky  shows  the  fit  of  the  data  to  the  model  [37] 

Sy  (f)  -Ky  /f(3  X  10-s  /Q4  ).  Ky  vahes  about  a  factor  of  500  for  Q-factors  between  10^  and 

3.8x10^  (resonator  frequencies  between  2.5  and  100  MHz). 

The  curve  labeled  K^  shows  the  fit  of  the  same  resonator  data  to  the  model  [38] 

S  (f)  -  K^  (3  X  10^0  /CK ).  K^  varies  about  a  factor  of  10  for  the  same  range  in  Q-factor.  Curves  Pe 
and  Pb  show  the  fit  of  the  same  resonator  data  to  the  model  Sy(f)  «  p/f  (Voi/Q^  ),  where  pe  is  for  SC 
and  AT  resonators  with  electrodes  plated  on  the  resonator  and  Pb  is  for  BVA-styie  AT  and  SC  resonators 
[43].  Volume  between  the  electrodes  (in  cm  3)  is  used  to  approximate  the  volume  of  quartz 
contributing  to  the  output  power.  The  B  factors  are  remarkably  constant  for  Q-factors  from  IQS  to  3.8 

X  10>. 

Figure  5  shows  the  dependence  of  Pe  on  Q-factor  for  a  number  of  electroded  resonators  of  the 
same  type  frcmi  a  single  manufacturer  for  3  resonator  types  as  measured  by  Norton  [44,51].  The  wide 
variation  in  p«  for  the  same  style  resonator  and  Q-factor  indicates  that  acoustic  loss  is  not  the  only 
mechanism  contrfouting  to  the  noise  level.  The  data  for  this  graph  was  taken  from  measurements  of 
(100  s)  ttfid  may  have  been  biased  high  by  random  walk  FM  noise  in  some  resonators. 

IV.  4.  Dfocussion 

The  analysis  of  the  most  stable  quartz  resonators  indicates  that  the  1/f  frequency  noise  level 
depends  on  volume  between  the  electrodes  and  unloaded  Q-factor  in  relatively  good  agreement  with  eq. 
(30)  considering  the  fact  that  the  estimation  of  <tii>  and  N  is  only  approximative.  The  nonadjustable 
parameters  p»  and  Pb  are  virtually  constant  versus  unloaded  Q-factor,  which  is  in  stark  contrast  to 
ttUng  parameters  K^  md  Ky.  It  is  not  surprising  to  us  that  p«  and  Pb  are  different  for  the  two  types  of 
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resonators,  sines  energy  trapping  and  electrode  stress  are  considerably  different.  Fig.  5  shows  that 
tttere  are  other  noise  processes  besides  acoustic  losses  that  affect  the  1/f  noise  level  in  some 
resonators. 

Although  we  have  analyzed  only  the  data  (br  a  few  resonators,  the  cor^ency  of  ^  and  Pb  over  a 
ta^or  of  40  in  Q-factor  and  resonator  frequency  and  the  general  agreement  for  the  magnitude  of  3 
between  theory  and  experiment  give  us  some  confidence  that  this  new  model  can  be  used  to  fxecSet  the 
best  performance  of  (Afferent  resonator  geometries. 

This  new  volume  model  predicts  that  a  resonator  having  smaller  electrodes  would  have  a  lower 
level  of  1/f  frequency  nofoe  than  another  one  with  the  same  frequency  and  Q-factor  but  larger  diameter 
eiectiodes.  The  decrease  in  electrode  area  would  increase  the  impedance  levels  and  degrade  the  wide¬ 
band  noise  somewhat.  For  most  resonators  the  wideband  noise  is  dominated  by  the  electronics  and  not 
the  resonator.  The  increase  in  series  resistance,  obtained  by  decreasing  the  eiectrode  area  by  a  factor 
of  4,  would  probably  be  tolerable  from  the  standpoint  of  wideband  noise  but  might  require  a  change  in 
loop  gain. 

BT  resonators  are  potentially  interesting  in  that  they  offer  a  Qvo  product  approximately  three 
times  higher  than  that  of  AT  and  SC  resonators.  BT  cuts  are  roughly  as  sensitive  to  temperature 
transients  as  AT  cuts.  Therefore  to  achieve  parts  in  frequency  stability  with  BT  cuts  would 
require  temperature  stabilities  of  order  10-*  K/s  or  100  times  better  than  is  required  for  SC  cut 
resonators  {51]. 

Based  on  these  early  observations  it  appears  that  the  level  of  1/f  frequency  noise  in  quartz  may 
yet  be  improved  to  the  low  10*^^  level  by  applying  one  or  more  of  the  following  techniques:  reducing  the 
eiectrode  area,  using  BVA  type  resonators,  going  to  lower  frequencies,  using  BT  cut  resonators,  it  must 
be  remembered  that  acceleration  induced  effects  become  more  dominant  as  the  stability  improves. 


V.  PRACTICAL  APPUCATION  TO  QUANTUM  1/f  NOISE  IN  INFRARED  DETECTORS 

Quantum  1/f  Noise  is  a  fundamental  aspect  of  quantum  mechanics,  representing 
universal  fluctuations  of  physical  process  rates  R  and  cross  sections  o  given  by  the  fractionai 
(or  relative)  spectral  density  S(f)  «  2aA/fN.  Therefore  it  is  present  in  the  process  rates 
generating  the  dark  current  observed  in  junction  photodetectors,  such  as  diffusion  (scattering 
cross  sections  fhictuate)  in  diffusion-limited  junctions,  and  recombination  in  the 
reoombfriation-limited  regime.  One  is  therefore  fompted  to  expect  similar  fluctuations  in  the 
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photogoneratton  of  etoctron-hoie  pairs.  However,  as  we  show  below,  the  oorresporxling  c^antum 
1/f  coeffidefit  is  zero,  precluding  tfie  existence  of  quantum  1/f  fluctuations  in  the  photo¬ 
generation  rate.  Here  N  is  the  number  of  carriers  used  to  define  or  measure  the  process  rate  or 
cross  section  considered. 

For  an  arbitrary  process  involving  a  total  of  n  incoming  and  outgoing  charged  particles, 
the  nonrelativistic  quantum  1/f  coefficient  is  given  [52]  by 

n 

2oA  -  (4a/3jic2)  XTtmjqiqj{vi-Vj)2.  (32) 

i:i-l 


where  the  summation  runs  over  the  charges  qi  and  velocities  vj  of  ail  incoming  (qi— 1)  and 
outgoing  (qi-l)  particles  (altogether  n  of  them)  in  the  process  whose  quantum  1/f  noise  we 
want  to  find,  and  a  is  Sommerfeid's  fine-structure  constant,  e2/Ko1/137.  In  a  photoelectric 
process  a  photon  (q«0)  is  absorbed,  and  a  pair  of  oppositely  charged  particles  is  generated 
(q-1)  with  velocities  vi  and  V2  which  are  either  zero,  or  quickly  decay  to  zero  in  a  time 
negligible  with  respect  to  the  reciprocal  frequency  at  which  we  calculate  the  quantum  1/f  noise. 
Thus  in  our  case  there  are  no  incoming  charged  particles,  and  n«0+2«2.  The  oA  coefficient  of  a 
photogeneration  process  is  therefore  zero, 

aAph  -  (1.1)+(2,2)+(1.2)+(2.1)  -  0+0+  (4o/3jtc2)(vi -V2)2  -  0.  (33) 

Ail  photogenerated  carriers  of  the  right  sign  are  collected  in  the  well  of  the  charge-coupled 
device,  although  tfwy  may  generate  quantum  1/f  voltage  fluctuations  on  their  way.  Since 
usually  only  the  number  of  carriers  collected  at  read-out  matters,  no  quantum  1/f  noise  wiil  be 
observed  in  a  photoelectric  CCO  as  long  as  ttie  dark  current  is  negligtole  with  respect  to  the 
photocurrent.  This  is  in  agreement  with  the  experiments  performed  by  Mooney  [53]  of  RADC- 
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H«isoofn  AFB.  The  same  considerations  apply  to  Metai-lnsulator-Semk»nductor  (MIS) 
photoctetectors  [54]. 

VI.  QUAivruM  1/r  mobility  fluctuations  in  homogeneous  semiconductor  samples 

VI.  1.  Introduction 

A  first  principles  calculation  of  quantum  1/f  cross  correlations  performed  [55]  for  the 
first  time  in  1987  by  one  of  us  [56]  has  yielded  a  slightly  different  result  compared  to  earlier 
expectations.  This  same  new  form  of  the  quantum  1/f  cross  correlations  was  derived  also  with 
a  (Afferent  method  by  Van  Vliet  in  1989.  It  differs  from  the  old  form  used  in  the  1985 
calculation  of  Kousik  et.  al.  by  a  correction  which  is  zero  when  the  momentum  changes  of  the  two 
current  carriers  involved  in  the  cross  correlation  are  identical,  but  increases  to  finite  values 
when  the  momentum  differences  caused  by  the  scattering  process  are  different.  The  correction 
is  proportional  to  the  squared  difference  of  the  two  momentum  changes.  We  have  repeated  all 
calculations  in  the  original  paper  by  Kousik  et.ai.  [56],  obtaining  both  for  impurity  scattering 
and  for  the  various  types  of  phonon  scattering  new  analytical  expressions  which  show  a 
considerable  increase  of  the  final  quantum  1/f  noise.  The  results  obtained  are  applicable  both  to 
(flrect  and  indirect  bandgap  semiconductors. 

We  have  performed  an  analytical  calculation  of  mobility  fluctuations  in  silicon  and 
gallium  arsenide,  using  the  new  quantum  1/f  cross-correlations  formula.  This  calculation  is  of 
major  importance  for  the  1/f  noise-related  optimization  of  the  two  types  of  materials,  and  of 
the  many  devices  constructed  with  them  for  military  and  civilian  applications  in  the  electronic 
and  opto-electronic  industry. 

The  new  cross-correiation  formula  gives  the  cross-spectral  density  which  describes  the 
way  in  which  simultaneous  quantum  1/f  scattering  rate  fluctuations  AW  observed  in  the 
direction  of  the  outgoing  scattered  wave-vector  K'  are  correlated  with  those  in  the  K”  direction, 
when  the  two  corresponding  incoming  current  carriers  have  the  wave  vectors  Ki  and  K2: 

SAW(Ki.K*;K2,K-:f) 

-  (2o/3xf)(IVm*C)2WKi.K*WK2,K-[(K*-Ki)2+(K--K2)2]8Ki,K2.  (34) 

The  form  conjectured  by  us  earlier  had  2(K'-Ki)(K”-K2)  in  place  of  the  rectangular  bracket. 
The  difforence  between  the  rectangular  bracket  and  2(K'-Ki)(K”-K2)  is  the  perfect  square 
[(K'-Ki)-(K”-K2)]2.  Therefore  we  expect  the  new  results  to  be  always  larger  than  the 
results  oMained  on  the  basis  of  the  previously  conjectured  form. 


22 


VI.  2.  Impurity  Scattering 


For  in^rity  scattering  of  electrons  in  solids,  fluctuations  Ax  of  the  collision  times  x 
will  cause  mobility  fluctuations 

Auband(t)  -  Ie/m*«v2»]ZKVK2Ax(t)nK.  (35) 

where  «v2»  is  both  the  average  over  all  states  of  wave-vectors  K,  with  occupation  numbers 
nK,  in  the  conduction  band,  and  the  thermal  equilibrium  average  of  the  quadratic  carrier 
velocities.  With  the  help  of  the  relation 

1/x(K)  -  (V/8ic3)J(1  -  COS0Vcose)WK,K*ri^K\  (36) 

the  mobility  fluctuations  are  reduced  to  fluctuations  of  the  elementary  scattering  rates  Wk,K'> 
governed  by  Eq.  (34).  Here  V  is  the  volume  of  the  nonnaiization  box  which  disappears  in  the 
final  result.  0  and  0'  respectively  the  angles  K  and  K*  form  with  the  direction  of  the  applied 
field.  One  finally  obtains  after  tedious  multiple  integrations 

^■2SA^(f)  -  (256KaKVhl2/3m*824e8Nj2](i/f)x,^Ki0lln(1+a2).a2/(l+a2)]- 
3l(2a2+a4)/(1+a2).2ln(1+a2)]F(EK)[lKVK2x(K)F(EK)]-2,  (37) 

where  a-2K/K,  K2«e2n(T)/ckBT,  n(T)  is  the  electron  concentration,  F(Ek)-oxp(Ef-Ek)  for 
non-degenerate  semiconductors,  Nj  the  concentration  of  impurities  of  charge  Ze  and  e  the 
dielectric  constant.  The  corresponding  partial  Hooge  parameter  for  impurity  scattering  is  thus 

oo 

ai  ■  [4V2jtaKn5Nc/3m*^^2(kgj)3/2c2j  ^^20-x 

0 

[In(bx+1)-bx/(bx+1)]-3((2bx+b2x2)/(bx+1)-2ln(bx+1)l 

PO 

{ Jdxx3e-xiln(bx+l)-bx/(bx+l)l-^}*2.  (38) 

0 

This  result  is  graphed  below  for  three  different  values  of  the  donor  concentration  Nd  and  is 
compared  with  old  results  obtained  by  simply  recaicuiating  the  old  analytical  expression  [56] 
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Alpha  impurity  New  and  Old 


Nd  -10E21  new 
Nd  -10E23  new 
Nd>10E24new 
Nd.10E21  old 
Nd-10E23old 
Nd  -10E24  old 


As  expected,  the  new  cross  correlation  formula  leads  to  higher  ai  values  than  the  previously 
conjectured  expression.  This  was  mentioned  in  connection  with  Eq.  (34)  above. 

Vi.  3.  Acoustic  Electron-Phonon  Scattering 

In  this  case  the  calculation  is  similar,  and  leads  to  the  result 

oo 

aac  -  [32*aNcm*C7h3/3c2kBT)^]{(1/R2)  Jdxx*-* 

1 

I(x-1)7/7  +  (R  +  1)(x-1)8/6  +  R(x-1)5/5] 

[(x-1)5/5  +  (R  +  1)(x-1)^/4  +  R(x-1)3/3]exp(-x2/4R) 

1 

+  Jdxx-*[(x  +  1)5/5-(x+1)6/6  +  (x-1)5/5  +  (x-1)6/6] 

0 

[(x  +  1)3/3  +  (x-1)^/4  +  (x-1)3/3-(x  +  1)4/4]exp(-x2/4R) 

+  Jdxx-<I(x+1)5/5-(x+1)6/6][(x-»-1)3/3-(x+1)^/4]exp(-x2/4R)}.  (39) 

1 

where  R-kBT/2m‘Ci2,  Ci  is  the  deformation  potential,  and  Nc  is  the  effective  density  of  states 
for  the  conduction  band. 
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Vt.  4.  Non-Poiar  Optical  Phonon  Scattering 


This  time  one  obtains 

PO 

on.o.ph  ■  [8irV2h(aoaNch2/3m*®^2c2j,j^j|  /dxx^^2 

0 

[(F+1)(x-1)1'20(x.1)  +  F(x  +  1)^'2j-4 

[(F+1)2(x-1)(2x-1)e(x-1)  +  F2(x+1)(2x+1)]exp(-ria)oX/kBT)} 

PO 

{ Jdxx3/2[(F4.^)(x-1)1/20(x.i)  +  F(x+1)^/2]-1exp(-riQ)oX/kBT)}*2,  (40) 

0 

where  F-[exp(n«i>o/i(BT)-1]*^.  and  (Dq  is  the  optical  phonon  frequency. 

VI.  5.  Polar  Optical  Phonon  Scattering 

Proceeding  as  in  Secs.  2  and  4.  we  obtain 

oo 

Op.o.ph  ■  [8icV2h<B|aNch2/3m*5/2c2{jj|]^  /dxx^ 

0 

(F2(x  +  1)1/2|n(2x1/2+2(x+1)1/2) 

+  (F+1  )2(x-1  )1/2|n(2x^/2+(x.i  )l/2)0(x-i  )]exp(-hfi)|x/kBT) 

((F+1)arcsinh(x-1)^^20(x-i)+Farcslnh(xi/2^j-4j_  (4.J  ^ 

Here  (di  is  the  longitudinal  phonon  frequency. 

Vi.  6.  Intervaliey  Scattering 

This  type  of  scattering,  present  in  indirect  bandgap  semiconductors,  transfers  electrons 
from  one  of  the  six  minima  (or  valleys)  of  the  conduction  band  energy  in  k-space  to  one  of  the 
other  five  minima.  Transitions  between  a  valley  and  the  nearest  valley,  which  is  along  the  same 
k-space  direction  in  the  next  copy  of  the  first  Briiloin  zone  in  the  periodic  zone  scheme,  are  of 
the  Umkiapp  type,  and  are  called  g-processes.  Trsmsitions  to  the  four  valleys  present  in  the 
same  zone  along  the  other  two  k-space  directions  are  called  f-processes.  Repeating  a  previous 
calculation  [56]  on  the  basis  of  the  new  aoss-correlation  formula  (34),  we  obtain  fo  g- 
processes 

Og  «  [8KV2h<0|joNcH2/3m'5/2c2<njj]^  J(jxx5/2 
[(F+1)(x-1)1'20(x.i)  +  F(x+1)1'2]-4 

[(F  +  1)2(x-1)(2x-1)0(x-1)  +  F2(x  +  1)(2x+1)]exp(-rKDijx/kBT)) 
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(42) 


{  Jdxx3/2{(F+i)(x.i)l/20(x-i)  +  F(x+1)^'2]-1exp(.n(OiiX/kBT)}-2, 

0 

where  hoaij  is  the  phonon  energy  corresponding  to  the  momentum  difference  required  by  the 
intervailey  transition.  For  the  corresponding  f-process  we  obtain  [57] 

Of  -  (ko/qo)2ag,  (43) 

where  ko/qo  >s  the  ratio  between  the  position  vector  of  a  conduction  band  energy  minimum  in  k 
space,  and  twice  the  distance  of  the  minimum  from  the  Briilouin  zone  boundary,  0.85/0.3  for 
siiicon.  There  are  three  g-type  aiphas  ogi,  ag2  and  ogs  (from  LA,  TA  and  LO  phonons 
repectiveiy)  and  three  f-type  contributions  afi,  af2  and  afs  (from  TA,  LA  and  TO  phonons). 
Their  vaiues  are  given  in  the  graph  beiow  and  are  a  few  times  larger  than  the  old  values. 

Alpha  for  Intervailey  Scattering  Processes 


Alpha  g1 
Alpha  g2 
Alpha  g3 
Alpha  f1 
Alpha  f2 
Alpha  f3 


The  various  quantum  1/f  contributions  derived  here  can  be  approximately  superposed  to 
yield  the  resultant  quantum  1/f  coefficient  according  to  the  rule 

OH  -  2i(|i/iii)2  oi  (44) 
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